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We discuss some problems of dissipative chaos for open quantum systems 
in the framework of semiclassical and quantum distributions. For this goal, 
we propose a driven nonlinear oscillator with time-dependent coefficients, i.e. 
with time-dependent Kerr-nonlinearity and time-modulated driving field. This 
model showing both regular and chaotic dynamics in the classical limit is re- 
alized in several experimental schemes. Quantum dissipative chaos is analyzed 
on the base of numerical method of quantum trajectories. Three quantities are 
studied: the Wigner function of oscillatory mode from the point of view of 
quantum-assemble theory and both semiclassical Poincare section and quan- 
tum Poincare section calculated on a single quantum trajectory. The compara- 
tively analysis of these distributions for various operational chaotic regimes of 
the models is performed, as well as scaling invariance in dissipative chaos and 
quantum interference effects assisted by chaos are discussed. 

Keywords: Quantum optics. Nonlinear dynamics 



1. Introduction 

Quantum nonlinear systems with chaotic classical counterparts have re- 
ceived much attention in the last two decades. This field of investigation 
is sometimes called quantum chaoi^^. The majority of studies of quantum 
chaos for isolated or so-called Hamiltonian systems, the classical counter- 
parts of which are chaotic, focus on static properties such as spectral statis- 
tics of energy levels and transition probabilities between eigenstates of the 
system. A variety of studies have also been carried out to understand the 
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features of time-dependent chaotic systems. In contrast to that very ht- 
tle work has been done to investigate quantum chaos for open nonhnear 
systems. The beginning of study of an open chaotic system can be dated 
back to the papertP where the authors have analyzed the kicked rotor and 
similar systems with discrete time interacting with a heat bath. Quite gen- 
erally, chaos in classical conservative and dissipative systems with noise has 
completely different properties, e.g., strange attractors can appear only in 
dissipative systems. For dissipative systems Poincare section has form of 
strange attractor in phase space while for Hamiltonian systems it has form 
of close contours with separatrices. 

In classical mechanics a standard characterization of chaos might be 
given in terms of the unpredictability of phase-space trajectories or Poincare 
section that consists of non-localized distributed points. However, the most 
important characteristic of classical chaotic systems-exponential divergence 
of trajectories, starting at arbitrarily close initial points in phase space 
does not have quantum counterpart because of the Hiesenberg uncertainty 
principle. The question has been posed of what constitutes the quantum 
mechanical equivalent of chaos. Many criteria have been suggested to de- 
fine chaos in quantum systems, varying in their emphasis and domain of 
applicatiorPJ As yet, there is no universally accepted definition of quan- 
tum chaos. 

In recent years much effort has been expended, both theoretically and 
experimentally, to explore the role of quantum fluctuations and noise in the 
order-to-chaos transition for open systems. It is obvious that the investiga- 
tions in this area are connected with the quantum-classical correspondence 
problem, in general, and with the environment induced decoherence and 
dissipation, in particular. Recently, this topic has been the focus of theoret- 
ical investigations. As a part of these studies, it has been recognized 131 that 
the decoherence has rather unique properties for systems classical analogs 
of which are chaotic. In particular, the formation of sub-Planck structure 
in phase space has been discussed for chaotic systerrP. The connection be- 
tween quantum and classical treatments of chaos was also realized by means 
of comparison between strange attractors on the classical Poincare section 
and the contour plots of the Wigner functionJISl. 

In this paper we investigate some problems of dissipative chaos for open 
quantum systems in the framework of quantum distributions. We use the 
traditional ensemble description of Markovian open systems, based on the 
master equation. This equation is presented in quantum trajectories in the 
framework of the quantum state diffusion approactP. Recently, it was shown 
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how quantum state diffusion can be used to model dissipative chaotic sys- 
tems on individual quantum trajectorieJ^ In contrast with these papers, 
here we show that it is possible to describe quantum chaos using also a 
statistical ensemble of trajectories, which is actually realized in nature. 

The requirement in realization of this study is to have a proper quantum 
model showing both regular and chaotic dynamics in the classical limit. We 
propose a driven nonlinear oscillator with time-dependent coefhcients for 
this goal. This model proposed to study the quantum chaos in the series of 
papers ^ allow us to examine challenging problems of quan- 

tum dissipative chaos, including the problem of the quantum counterpart 
of a strange attractor. 

The other problem of our interest in this paper relates to the quantum 
effects in systems with chaotic dynamics. Particularly, in a recent paper, 
sub-Poissonian statistics of oscillatory excitations numbers was established 
for chaotic dynamics of nonlinear oscillatoJ^. It was shown that quantum- 
interference phenomena can be realized for the dissipative no nlinear sys- 
tems exhibiting hysteresis-cycle behavior and quantum chaoJ^ The 
study of these phenomena provides a fundamental understanding of quan- 
tum fluctuations in quantum chaos and opens a way for new experimental 
studies of the quantum dissipative chaos in the field of quantum optics. 

The outline of this paper is as follows. In Sec. II we describe both the 
model and the method of calculations. In Sec. Ill we analyze correspondence 
between Poincare sections and Wigner functions for the chaotic dynamics. 
In Sec. IV we analyze scaling invariance for quantum system and discuss 
quantum chaos for the regimes in which the classical chaos is lost. In Sec.V 
we present results on quantum interference phenomena for chaotic dynam- 
ics. We summarize our results in Sec. VI. 



2. Driven nonlinear oscillator as an open quantum system 

In this section the systems and the methods of calculations are presented. 
We treat the Duffing oscillator as an open quantum system and assume 
that its time evolution is described by Markovian dynamics in terms of the 
Lindblad master equation for the reduced density matrix p. In the interac- 
tion picture that corresponds to the transformation 

where a"*" and a are the Bose annihilation and creation operators of the os- 
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cillatory mode and lo is the driving frequency, this equation reads as 

^ (upL+ - ^L+Up - ^pL+L^ . (1) 

The Hamiltonians are 

Hq = hAa'^a, 

Hm = ftx(0(«+a)' + Kf{t)a+ + fit)* a), (2) 

where x(^) ^-nd f{t), which may or may not depend on time, represent, 
respectively, the strength of the nonhnearity and amplitude of the force, wg 
is the resonant frequency, A = luq — uj is the detuning. The dissipative and 
decoherence effects, losses, and thermal noise are included in the last part 
of the master equation, where Li are the Lindblad operators: 

Li = ^/{N + l)ja, L2 - ^/N^a+, (3) 

7 is the spontaneous decay rate of the dissipation process and N denotes 
the mean number of quanta of a heat bath. Here we focus on the pure 
quantum effects and assume = 0. 

This model seems experimentally feasible and can be realized in several 
experimental schemes. In fact, a single mode e.m. field is well described in 
terms of an anharmonic oscillator (AHO), and the nonlinear medium could 
be an optical fiber or a x(3) crystal, placed in a cavity. The anharmonicity 
of mode dynamics comes from the self-phase modulation due to the photon- 
photon interaction in the x(3) medium. In this case, it is possible to realize 
time modulation of the strength of the nonhnearity by using a media with 
periodic variation of the x(3) susceptibility. 

On the other side, the Hamiltonian described by Eq. ([2|) describes a 
single nanomechanical resonator with a'^ and a raising and lowering oper- 
ators related to the position and momentum operators of a mode quantum 
motion 

X = \ (a + a^), p — ~i\/2hm(jJo(a — a^), 

V zmujo 

where m is the effective mass of the nanomechanical resonator, ujq is the lin- 
ear resonator frequency and x proportional to the Duffing nonhnearity. One 
of the va riants of nano-oscillators is based on a double-clamped platinum 
beam for which the nonhnearity parameter equals to x = h/4:\/3Qma^, 
where Oc is the critical amplitude at which the resonance amplitude has 
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an infinite slope as a function of the driving frequency, Q is the mechan- 
ical quality factor of the resonator. In this case, the giant nonlinearity 
X = 3.4 • 10~'*s~^ was realized. Note, that details of this resonator, includ- 
ing the expression for the parameter Oc, are presented i On decreasing 
nanomechanical resonator mass, its reso nance frequency increases, exceed- 
ing 1 GHz in recent experimentJ^^ ^\ \i jg possible to reach a quantum 
regime for such frequencies, i.e., to cool down the temperatures for which 
thermal energy will be comparable to the energy of oscillatory quanta. 
The recent investigations in this direction are devoted to classical to quan- 
tum transition of a driven nanomechanical oscillatoj ^, g eneration of Fock 
stateJl^, nonlinear dynamics, and stochastic resonance^^. Note, that in the 
last few years there has been rapid progress in the construction and manipu- 
lation of such nanomechanical oscillators with giant x(3)-Kerr nonlinearity. 
The nanomechanical resonator with a significant fourth-order nonlinearity 
in the elastic potential energy has been experimentally demonstrated ^^II \i 

has also been shown that this system is dynamically equivalent to the Duff- 
fool 

ing oscillator with varied driving force This scheme is widely employed 
for a large variety of applications as well as the other schemes of micro- 
and nanomechanical oscillators, more commonly as sensors or actuators in 
integrated electrical, optical, and optoelectrical systems-*^ . 

Cyclotron oscillations of a single electron in a Penning trap with a mag- 
netic field are another realization of the quantum version of the Duffing 
oscillatoJ^^^.In this case the anharmonicity comes from the nonlin- 
ear effect that is caused by the relativistic motion of an electron in a trap, 
while the dissipation effects arise from the spontaneous emission of the syn- 
chrotron radiation and thermal fluctuations of the cyclotron motion. Note 
that a one-electron oscillator allows one to achieve a relatively strong cubic 
nonlinearity, xll ^ 1- 

In recent years the study of quantum dynamics of oscillators with time- 
dependent parameters has been focus of considerable attention. This in- 
terest is justified by many applications in different contexts. Particulary, 
one application concerns to the center of mass motion of a laser cooled 
and trapped ion in a Paul trajj^. The quantum dynamics of an AHO 
with time dependent modulation of its frequency and nonlinearity param- 
eters has been inve stigated in applications to macroscopic superposition of 



It is well assessed that in the case of unitary dynamics, without any 
losses, an anharmonic oscillator leads to sub-Poissonian statistics of oscil- 
latory excitation number, quadratic squeezing and superposition of macro- 
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scopically distinguishable coherent states. For dissipative dynamics the im- 
portant parameter responsible for production of nonclassical states via x(3) 
materials is the ratio between nonlinearity and damping. Therefore, the 
practical realization of such quantum effects requires a high nonlinearity 
with respect to dissipation. In this direction the largest nonlinear interac- 
tion was proposed in many papers, particularly, in terms of electromag- 
netically induced transparencjE^ and by using the Purcell effeclPSl, and 
in cavity The significant nonlinearity has also been observed for 

nanomechanical resonatorJSI. These methods can lead to x(3) nonlinearity 
of several orders of magnitude higher than natural optical self-Kerr in- 
teractions. Note, that high x(3) nonlinear oscillators generate also a lot of 

lOOl 

interest recently due to their applications in areas of quantum computing^. 

In the case of nonlinear dissipative x(3) interaction stimulated by co- 
herent driving force, the time evolution cannot be solved analytically for 
arbitrary evolution times and suitable numerical methods have to be used. 
Nevertheless, with dissipation included a driven AHO model has been solved 
exactly in the steady-state regime in terms of the Fokker-Planck equation 
in the complex P representation- -'. Analogous solution has been obta ined 
for a combined driven parametric oscillator with Kerr nonlinearitj031. The 
Wigner functions for both these models have been obtained using these 
solutionJ^Sl [35| 

The investigation of quantum dynamics of a driven dissipative nonlinear 
oscillator for non-stationary cases is much more complicated and only a few 
papers have been done in this field up to now. More recently, the quantum 
version of dissipative AHO or the Duffing oscillator with time-modulated 
driving force has been studied in the series of the papers ' in the context of 
a stochastic resonanc qu antum-t o-classical transition and investigation 
of quantum dissipative chaoJI^ 

For the constant parameters x(t) — x and f{t) = f the equations ([1]) 
and ([2]) describe the model of a driven dissipative AHO that was introduced 
long ag o in quantum optics to describe bistability due to a Kerr nonlinear 
mediunJSfiJ. For the case of time-dependent parameters x{t) and f{t) the dy- 
namics of the AHO exhibits a rich phase-space structure, including regimes 
of regular, bistable and chaotic motion. We perform our calculations for reg- 
ular and chaotic regimes concerning two models of time-modulated AHO 
corresponding to two physical situations: (i) x = x{t) = Xo + Xisin{i}t) 
and f{t) = const = /; (ii) / = f{t) = fo + fi exp{St) and x{t) = const = x 
with S <^ uj and fl <^ uj are the modulation frequencies. 

We analyze the master equation numerically using quantum state dif- 
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fusion method According to this method, the reduced density op- 

erator is calculated as an ensemble mean 

1 ^ 

p(i) = M(|^5W)(^5(t)|)= lim J2\i:^it)){Mt)\ (4) 

over the stochastic pure states \ip^{t)) describing evolution along a quan- 
tum trajectory. The stochastic equation for the state |V'c(i)) involves both 
Hamiltonian described by Eq. ^ and the Linblad operators described by 
Eq. and reads as: 

< U X L+ >)\^{)dt + J2 < L^ >)|*?)dC, (5) 

1=1,2 

where is the generated complex gaussian noise that satisfies the following 
conditions: 

M(dCO - 0, M{di^dij) = 0, MididCj ) = ^^odt. (6) 

We calculate the density operator using an expansion of the state vector 
\'ip() in a truncated basis of Fock's number states of a harmonic oscillator 

n 

3. Wigner function and Poincare section 

In this section we shortly examine correspondence between the Poincare 
section and the Wigner function of the oscillatory mode considering the 
model (ii) with time-dependent driving amplitude In semiclassical 

treatment a chaotic operational regime is analyzed on phase space of di- 
mensionless position and momentum x = Re{a) and y = Im{a), where 
a —< a> is the oscillatory complex amplitude (see, Eq. ifTl]) '). Choosing xq 
and t/o s-s an arbitrary initial phase-space point of the system at the time 
tQ, we define a constant phase map in the plane by the sequence of points 
at tn = to + (2TT/S)n, where n — 0,1, 2.... This means that for any t ~ tn 
the system is at one of the points of the Poincare section. 

The analysis show that for time scales exceeding the damping time, 
t 3> 7~^, the asymptotic dynamics of the system is regular in the limits 
of small and large values of the modulation frequency, i.e., 5 ^ ^, 5 "f 
and small or large amplitude of driving field /i ^ /o, /i ^ /o- Fig. [TJc) 
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shows the results of numerical calculations of the classical maps, for the 
parameters chosen in the range of chaos, i.e. /o ~ /i and 5 > 7. As we 
see, the figure clearly indicates the classical strange attractor with fractal 
structure that is typical for a chaotic dynamics. 

It should be mentioned, that the Wigner function is one of the quantities 
that allows to observe chaos in quantum theory. The nonstationary Wigner 
function is written as for the state vector li^it))- 

Ty(a,i) = ^exp(-2|an J f3 f3\m) {mW) exp (-2(/3a* - /3*a)) 

(8) 

We apply the QSD to determine Wigner functions for the quantum 
states of a driven anharmonic oscillator during time evolution. For this, 
we use the well-known expression for the Wigner function in terms of the 
matrix elements pnm = {n\p\m) of the density operator in the Fock state 
representation: 

Wir,e) =yPnm{t)Wranir,0). (9) 



Here: (r, 0) are the polar coordinates in the complex phase-space plane, 
X = rcosO, y — rsin6', while the coefficients WmnirjO) are the Fourier 
transform of matrix elements of the Wigner characteristic function: 



f (-1)" y ^ exp i{m - n)0(2r)™-" 
exp(-2r2)i™-"(4r2), m > n 

I (-1)" V^expi(m ~ n)^^(2r)"-" '' 
exp(-2r2)L;57™(4r2), n > m. 



(10) 



It is remarkable that there is a correspondence between contour plots of 
the Wigner function and the Poincare section. This point is illustrated on 
the Fig.[Tl where the Wigner function Fig. (iKa), its contour plot, Fig.[TKb) 
and Poincare section, Fig.[TKc) are presented for the same parameters. Note, 
that the Wigner function is a quasidistribution in phase space averaging 
an ensemble of quantum trajectories obtained for a define time moment 
while Poincare section is the distribution for time intervals: it is constructed 
by fixing points in phase space at a sequence of periodic moments. What 
we can conclude from their correspondence is the fact, that the quantum 
distribution in phase space corresponds on the form to the semiclassical 
distribution but for big numbers of time intervals. As it is seen, the Poincare 
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Fig. 1. The Wigner function (a) and its contour plot (b) averaged over 3000 trajectories. 
The Poincare section (c) (20000 points plotted at times of the constant phase for the 
same parameters). In both cases, the dimensionless parameters are in the range of chaos, 
i.e. xh = 0.1, A/7 = -15, /o/7 = hh = 27, 5/7 = 5. 

section has fine fractal structure while Wigner function contour plot has not. 
This is due to the Hiesenberg uncertainty relations which prevent sub-Plank 
structures in phase space. 



4. Scaling invariance in dissipative chaos 

In this section the scaling invariance for the case of chaotic dynamics is 
considered on the base of anharmonic oscillator with time-dependent driven 
force. In the classical limit the system is described by the following equation 
of motion for the dimensionless amplitude: 

^ = -^a - z(A + x(l + 2|an)a - z(/o + cxp(-<5t)). (11) 

This equation is invariant for the scaling transformation of complex am- 
plitude a — \a if the other parameters transforms like: A — )• A = 
A + x(l - 1/A2), X ^ x' - x/A^ / ^ ./' = A/, 7 -> 7' = 7. This 
scaling property of the classical equation for the chaotic dynamics leads to 
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the symmetry of strange attractors: they have the same form in the phase 
space and differ from each other only in scale. This fact is demonstrated in 
the Fig. m 



(a) 



(b) 



lm(a) 0- 



lm(a) 




Re(a) 



-6 -4 -2 2 4 6 
Re(a) 



Fig. 2. The Poincare section for the following parameters: A/7 = —15, x/l = 2, 
/0/7 = 5.8, /1/7 = 4.9, 5/7 = 2, [2l;(a)); for the scaled A = 2 parameters [21 (b)). The 
Poincare sections are generated for time evaluation 5000 force periods. 



Strictly speaking the quantum system does not obey the same scaling 
invariance as classical one. An analysis of scaling invariance from the point 
of view of quantum-statistical theory has been performed!^. It was shown 
that such parameter scaling occurs for wider ranges of the parameters, but 
for not large values of the parameter x/7i where system displays strong 
quantum properties. 

Now we use the scaling arguments considering stochastic dynamic of a 
single trajectory in this way: to investigate chaos we use Poincare section 
based on evolution of a single quantum trajectory, Eql5l The Poincare sec- 
tion is obtained by recording {Re(a^), Im{a^)) at time intervals of 27r/5, 
where ^ indicates the stochastic variable and is obtained from Eq. [5l 
Our goal is to analyze the scaling invariance on the base of this quantity. 
On the other side we present results that relies on quantum-classical corre- 
spondence. 

In this direction, we consider the specific parameters for which DAO is 
in the vicinity of chaotic behavior determined classically. That means, if the 
parameters are slightly tuned in this range the transition from chaotic to 
regular dynamics might be realized. As we have realized above the dynamics 
of the system is chaotic in the range of parameters /o — /i and 5 > 7. Par- 
ticularly, it is demonstrated on the base of the semiclassical Poincare section 
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the system exhibits chaotic dynamics for /0/7 = /1/7 = 5.8, A/7 = —15 
and 5/7 = 2 in classical treatment. As analysis shows for the considered pa- 
rameters the system dynamics continues to be chaotic till /1/7 = 4.9(Fig[2), 
while chaotic dynamics becomes regular for /1/7 — 4.8. Thus, we examine 
Poincare section in quantum trajectory started from the regular regime, for 
the parameter /1/7 = 4.8. The results of calculations are presented in Fig. [3] 
for the parameters: A/7 = -15, x/7 = 2, /0/7 = 5.8, S/j = 2, /1/7 = 4.8. 
As we see, for the scaled parameter A = 1 the Poincare section is distributed 
stochastically but approximately exhibits localization in two ranges (Fig. 
[3l^a)); its shape does not coincide with Poincare section (Figl2][c)). For the 
case of A = 2, points of the Poincare section on a trajectory are distributed 
chaotically, however, unlike to the previous case A = 1, here an correspon- 
dence between both shapes take place. This likeness become more obvious 
for the case of A = 3 (Fig.[He)). 

It is important to note that in the semi-classical treatment all cases 
A = 1,2,3 correspond to regular regime as it is expected due to scaling 
invariance. Another situation is realized in the quantum treatment. For A = 
1 the system seems to be in regular regime. This statement is also confirmed 
by calculation of the Wigner function. The calculations shows that for such 
parameters the Wigner function has two-peak structure indicating that 
the regime of bistability is realized. If we increase the scaling parameter 
the resulted Poincare section's shapes in a quantum trajectory in their 
forms are coincided with Poincare sections in semiclassical treatment. So 
we found the parameters for which the system in classical treatment has 
regular dynamic while in quantum treatment its dynamic is chaotic. 

The results on mean excitation numbers are presented on Fig. [Sjb, d, 
f). As we see, for A = 1, the oscillatory excitation numbers varies from 1 
to 4 and thus the system is in deep quantum regime. The level of quantum 
noise is comparatively sufficient and chaos cannot exhibit itself on Poincare 
section. But for scaled parameters A = 2, 3 the excitation numbers increase 
and varies up to 40 (Fig.[3Id,f)). For these cases the ranges of variation are 
much enough to exhibit a chaotic type structures. 

Note, the usefulness of scaling procedure. It allows us to analyze clas- 
sically analogy regimes for various groups of the parameters. In particular, 
increasing scaling parameter of the system we can investigate quantum- 
classical correspondence. 
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Fig. 3. The Poincare section calculated on the base of a single trajectory for the follow- 
ing parameters: A/7 = -15, xH = 2, fo/f = 5.8, /1/7 = 4.8, S/-y = 2, (Fig. [3] (a)); the 
excitation number for the same parameters (Fig. |3jb)). Figs.|3jc) and[3ld) corresponds 
to scaled A = 2 parameters. Figs.lSfe) andlSff) corresponds to scaled A = 3 parameters. 
The Poincare sections are generated for time evaluation 5000 force periods. 
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5. Quantum interference assisted by chaos 

It is well assessed that in the case of unitary dynamics, without any losses, 
an anharmonic oscillator leads to sub-Poissonian statistics of oscillatory ex- 
citation number, quadratic squeezing, and superposition of macroscopically 
distinguishable coherent states. The dissipation and decoherence lead usu- 
ally to losses of these effects. In this section, we have numerically studied 
the phenomena at the overlap of chaos, dissipation, and quantum effects 
for the time-dependent nonlinear model. It was recently showrfUl 
that physical systems based on this model have a potential for genera- 
tion of high-degree sub-Poissonian light as well as for the observation of 
quantum-statistical effects and quantum interference that accompanied by 
chaotic dynamics. Here we concentrate on studies of quantum interference. 
We have pointed out that the time modulation of the oscillatory parame- 
ters, which are the strength of third-order nonlinearity or the amplitude of 
the driving force, leads to formation of the quantum-interference patterns 
in phase space in over transient regimes, for the definite time intervals ex- 
ceeding the transient dissipation time. 

It is well known that the phase-space Wigner distribution function can 
simply visualize nonclassical effects including quantum-interference. For ex- 
ample, a signature of quantum interference is exhibited in the Wigner func- 
tion by non-positive values. In this section the numerical results of the 
nonstationary Wigner functions in chaotic regimes of AHO are presented 
and discussed. 

It should be noted that the most of investigations of the quantum distri- 
butions of oscillatory states, including also modes of radiations, have been 
made for the steady-state situations. The simplicity of Kerr nonlinearity 
allows to determine the Wigner function of the quantum state under time 
evolution due to interaction. In this sense, we note the main peculiarity of 
our paper in comparison with above noted important inputs. In this paper, 
we calculate the Wigner functions in an over transient regime, t 3> 7~^, 
of the dissipative dynamics, however, we consider time-dependent effects 
which appear due to the time-modulation of the oscillatory parameters. 

Below we investigate the Wigner functions for time-modulated nonlin- 
earity. Note that the dynamics of the system with time-modulated nonlinear 
strength AHO is chaotic in the ranges: <5 > 7 and xo — Xi ^^'^ for neg- 
ative detuning. As shows analysis, controlling transition from the regular 
to chaotic dynamics can be realized through the intermediate ranges of 
bistability by varying the strength xi of the modulation processes in the 
ranges from xi Xo to Xi ^ Xo- The results of the ensemble averaged 
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-8 -4 4 8 



Fig. 4. The mean excitation number (a); the Poincare section (b) (si 20000 points) for 
the dimensionless complex ampHtude a, plotted at times of the constant phase St^ = 27rfc 
(k = 0,1,2...), when the maximal interference pattern on the Wigner function (c) is 
realized, for the case of time-modulated nonlinearity x(i)- The parameters are in the 
range of chaos: A/7 = 5,x(t)/7 = 0.2(1 -f 0.75 sin(m)), n/7 = 3J/-y = 10. 

numerical calculation of the mean excitation number, the Poincar section 
and the Wigner function are shown on the Figs. |4] respectively. The mean 
excitation number of the driven AHO versus dimensionless time is depicted 
in Fig. [4] (a). In contrast to the semiclassical result its quantum ensemble 
counterpart (see, Fig. [D^a)) has clear regular periodic behavior for time 
intervals exceeding the characteristic dissipation time, due to ensemble av- 
eraging. The Fig. IH^b) clearly indicates the classical strong attractors with 
fractal structure that are typical for a chaotic Poincare section. Thus, the 
Wigner function (Fig. |4l[a)) reflects the chaotic dynamics, its contour plots 
in the (x, y) plane are similar to the Poincare section. However, the Wigner 
functions have regions of negative values for the definite time intervals. The 
example depicted on Fig.[4l^a) corresponds to time intervals ^tk = 6+ ^f^7 
(fc = 0, 1, 2...), for which the mean excitation number reaches a macroscopic 
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level, i.e., n = 52. The interference pattern is destroyed as the time mod- 
ulation is decreased. Indeed, it is showrPS EH that for oscillatory mode of 
driven anharmonic oscillator including dissipation the Wigner function is 
positive in all phase space. 

6. Conclusion 

In summary, we have discussed dissipative chaos answering what is the 
counterpart of the semiclassical Poincare section in quantum treatment. We 
have presented a type of nonstationary systems showing chaotic dynamics 
and intrinsically quantum properties which are modelled by a driven dissi- 
pative anharmonic oscillator with time-dependent parameters. The connec- 
tion between quantum and classical treatments of chaos has been realized 
by means of a comparison between strange attractors on the semiclassical 
Poincare section, the shapes of the Poincare section on a single quantum 
trajectory and the contour plots of the Wigner functions. We have anal- 
ysis the borders of validity of scaling invariance for quantum dissipative 
chaos and we have demonstrated realization of long-lived quantum inter- 
ference assisted by chaotic dynamics for over transient regime and for the 
macroscopic level of oscillatory excitation numbers. 
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